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1.0  INTRODUCTION  AND  SUMMARY 


Coupling  between  the  combustion  process  and  the  acoustics  of  the  combus¬ 
tion  chamber  are  important  factors  in  determining  combustion  stability  of  a 
solid  propellant  rocket.  This  coupling  results  because  the  combustion  process 
reacts  to  both  the  local  acoustic  pressure  and  the  local  acoustic  velocity. 
Because  of  the  complexity  of  both  processes,  they  cannot  be  totally  character¬ 
ized  analytically;  therefore,  laboratory  test  data  are  needed  in  making  analyt 
ical  combustion  stability  predictions.  One  attractive  test  method  conceived 
by  CSD  is  the  rotating  valve  apparatus.  CSD  conducted  two  programs  under 
AFRPL  contract  to  develop  and  demonstrate  the  rotating  valve  method  of  mea¬ 
suring  the  pressure  coupled  combustion  response  of  solid  propellants.  The 
results  show  agreement  with  T-burner  measurements  when  the  T-burner  vent  term 
is  taken  to  be  zero.  In  addition,  reproducible  operation  of  the  apparatus  has 
been  demonstrated  at  pressures  up  to  1,500  psi  with  propellants  containing  as 
much  as  18J  aluminum. 

The  rotating  valve  apparatus  also  has  the  potential  for  measuring  the 
velocity  coupled  response  function  of  solid  propellants.  Velocity  oscilla¬ 
tions  of  controlled  frequency  and  amplitude  can  be  generated  by  simultaneously 
operating  a  rotating  valve  at  each  end  of  the  motor,  180°  out  of  phase.  In 
this  configuration,  velocity  coupling  dominates,  and  the  effects  of  other 
processes,  such  as  pressure  coupling,  are  minimized.  With  this  modification, 
the  rotating  valve  method  offers  the  potential  for  experimentally  and  quanti¬ 
tatively  investigating  many  characteristics  of  velocity  coupling  which  have 
been  postulated  by  purely  analytical  arguments.  The  nature  of  these  character 
istics  determines  the  manner  in  which  velocity  coupling  is  incorporated  into 
the  overall  combustion  stability  analysis  of  a  solid  propellant  rocket  motor; 
thus  experimental  evaluation  of  velocity  coupling  characteristics  is  essen 
tial. 


Under  AFOSR  contract  No.F49620-77-C-0048,  CSD  has  investigated  the  dual 
valve  approach  for  measuring  velocity  coupling  characteristics.  This  investi¬ 
gation  included: 


1-1 


A.  Studies  to  analyze  the  transient  ballistics  and  provide  a  method  for 
data  interpretation  and  analysis 

B.  Parametric  studies  to  evaluate  momentum  effects  and  determine  the  fre¬ 
quency  range  for  which  valid  data  can  be  obtained 

C.  An  error  assessment  of  the  methods  used  to  reduce  experimental  data 

D.  The  design  and  fabrication  of  experimental  apparatus 

E.  Cold-flow  studies  to  evaluate  the  performance  of  the  apparatus  and  to 
develop  methods  for  calibrating  the  apparatus  for  data  reduction 
purposes 

F.  Combustion  tests  to  evaluate  the  limits  of  the  apparatus  performance, 
to  determine  the  ability  to  distinguish  differences  in  velocity  re¬ 
sponse  functions,  to  study  the  effects  of  internal  flow  patterns  on 
the  derived  velocity  response,  and  to  establish  the  applicability  of 
the  data  analysis  procedures. 

Substantial  progress  was  made  towards  developing  the  dual  rotating  valve. 
Analytical  studies  developed  two  models  for  the  burner,  one  based  on  linear 
velocity  coupling  and  the  other  based  on  nonlinear  (or  amplitude  dependent) 
velocity  coupling.  Using  one  dimensional  gas  dynamics,  these  models  also 

include  the  effects  of  pressure  coupling,  flow  turning,  and  particle  damping. 

These  models  predict  the  behavior  of  the  burner  to  driving  produced  by  two 
valves  operating  180°  out  of  phase,  two  valves  operating  in  phase,  and  a  sin¬ 
gle  valve  at  either  end  of  the  motor.  Comparisons  show  these  models  predict 
the  proper  behavior  at  low  frequency,  the  proper  system  damping  at  the  natural 

frequencies,  and  agree  with  the  model  for  the  modulated  throat  motor. 

These  models  were  then  used  to  develop  methods  for  deriving  velocity 
response  functions  from  experimental  data.  Reasonable  accuracy  was  demon¬ 
strated  for  the  linear  velocity  coupling  model.  For  the  nonlinear  model, 
response  functions  can  be  derived  if  the  threshold  speed  is  known.  Separating 
the  threshold  from  the  response,  however,  is  extremelv  complex.  Progress  was 
made,  but  development  of  a  useful  data  reduction  procedure  was  not  completed. 

Cold-flow  tests  demonstrated  satisfactory  apparatus  performance.  Methods 
for  calibrating  the  mechanical  alignment  were  also  verified  in  these  tests. 
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Combustion  tests  demonstrated  satisfactory  apparatus  performance  with 
both  aluminized  and  nonaluminized  propellants.  The  data  show  reasonable 
reproducibility,  significant  differences  with  formulation  changes,  and  an 
effect  of  frequency.  Frequency  spectra  show  significant  harmonic  content  which 
appears  to  be  generated  by  the  combustion  process.  Comparisons  with  analyt¬ 
ical  predictions  suggest  the  nonlinear  velocity  coupling  model  is  qualitatively 
consistent  with  the  data. 

These  results  demonstrate  the  dual  rotating  valve  technique  warrants  fur¬ 
ther  study  as  a  method  for  investigating  velocity  response  of  solid  propel¬ 
lants.  This  is  especially  true  when  one  considers  the  early  stage  of  develop¬ 
ment  of  the  apparatus,  and  the  data  reduction  process,  as  well  as  the  lack  of 
sound  mechanistic  understanding  of  the  velocity  coupling  process.  These 
studies  should  evaluate  the  two  combustion  models  for  interpreting  velocity 
coupling  effects,  should  investigate  the  basic  mechanism  in  greater  detail, 
and  should  explore  the  correlation  between  laboratory  and  motor  results. 


2.0  TECHNICAL  ACCOMPLISHMENTS 


2.1  DESCRIPTION  OF  THE  DUAL  ROTATING  VALVE  METHOD  FOR  VELOCITY  RESPONSE 

The  basic  approach  in  the  dual  rotating  valve  method  is  to  mount  two 
rotating  valves  on  the  combustion  chamber,  one  at  each  end,  as  shown  in  figure 
2-1.  Conventional  nozzles  are  also  mounted  at  each  end  of  the  test  motor  to 
control  the  steady  state  pressure  and  the  location  along  the  propellant  grain 
at  which  the  mean  velocity  is  zero.  Thus,  the  instantaneous  nozzle  area  at 
each  end  of  the  motor  may  be  represented  by  the  sum  of  a  steady  state  compo¬ 
nent  and  an  oscillating  component  as  indicated  in  figure  2-1.  If  the  two 
valves  are  in  phase,  the  oscillating  components  add,  as  shown  in  figure  2-2, 
producing  an  oscillating  component  in  the  total  area.  The  resulting  behavior 
is  identical  to  the  single  rotating  valve  configuration  which  is  dominated  by 
pressure  coupling  effects.  If  the  valves  are  180°  out  of  phase,  the  area 
oscillation  produced  by  one  valve  exactly  cancels  that  produced  by  the  other 
valve  and  there  is  nc  net  area  oscillation  to  provide  pressure  coupling.  Sig¬ 
nificant  velocity  oscillations  are  produced  in  the  test  motor,  however, 
because  the  venting  of  combustion  gas  alternates  between  the  two  ends  of  the 
motor.  These  velocity  oscillations  couple  with  the  combustion  to  produce 
burning  rate  oscillations.  With  a  constant  net  vent  area,  these  burning  rate 
oscillations  produce  measurable  pressure  oscillations.  These  pressure  oscilla¬ 
tions  can  be  related  to  the  burning  rate  response  by  using  a  ballistic  anal¬ 
ysis  of  the  burner.  Thus,  by  operating  the  two  valves  180°  out  of  phase, 
velocity  coupling  can,  in  principle,  be  studied  under  conditions  which  maximize 
velocity  coupling  contributions  and  minimize  pressure  coupled  contributions. 

Eisel1  and  co-workers  suggested  a  similar  approach  to  generating 
controlled  velocity  oscillations  several  years  ago.  Their  objective  was  to 
observe  the  luminosity  response  of  aluminized  propellants  to  velocity 
oscillations.  Their  experiment  incorporated  two  identical  choppers  for 
modulating  the  exhaust  flow.  They  did  not,  however,  include  additional 
conventional  nozzles  for  independently  controlling  the  mean  flow  environment 
nor  did  they  attempt  to  measure  and  interpret  the  pressure  oscillations 
quantitatively. 
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Figure  2-1.  Valve  Driving 
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Figure  2-2.  Valve  Layout  for  Velocity  Response  Measurement 
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The  dual  rotating  valve  approach  to  measuring  the  velocity  coupled 
response  offers  three  potentially  significant  advantages  over  the  T-burner 
method.  In  the  T-burnerf  the  test  samples  are  located  where  both  velocity 
oscillations  and  pressure  oscillations  produce  burning  rate  oscillations.  In 
fact,  the  maximum  driving  for  the  velocity  oscillations  is  at  best  only  50$  of 
the  total  driving.  Since  the  T-burner  measures  the  total  burning  rate 
response,  the  pressure  coupled  response  must  be  subtracted  from  the  total 
response  to  derive  the  velocity  response.  This  requires  independent  measure¬ 
ment  of  the  pressure  response.  More  importantly,  the  uncertainties  in  the 
total  measured  response  and  pressure  response  add  together  to  increase  the 
uncertainty  in  the  velocity  response.  The  dual  rotating  valve  method,  however, 
maximizes  the  velocity  driving  and  minimizes  the  pressure  driving,  which  should 
improve  the  accuracy  of  the  velocity  response  results  significantly. 

A  second  potential  advantage  involves  the  cost  of  testing.  The  T-burner 
method  requires  10  to  20  tests  at  each  condition  of  pressure  and  frequency  to 
derive  a  velocity  response  function  value.  In  the  rotating  valve  method,  only 
one  test  is  required  to  obtain  a  response  function  value  at  each  condition. 
There  is  also  a  corresponding  reduction  in  the  quantity  of  propellant  required, 
which  is  an  important  consideration  when  several  alternative  propellant  formu¬ 
lations  are  being  considered  for  a  particular  motor  application. 

A  third  potential  advantage  is  the  capability  for  varying  the  mean  flow 
environment.  Current  models  for  velocity  coupling  postulate  that  the  mean 
flow  environment  may  have  a  significant  effect  on  the  velocity  response  func¬ 
tion.  In  the  dual  rotating  valve,  shifting  mean  nozzle  area  from  one  end  of 
the  burner  to  the  other  end  changes  the  axial  velocity  profile  in  the  test 
motor.  Thus,  the  response  can,  in  principle,  be  measured  under  different  mean 
flow  environments.  The  T-burner  does  not  permit  this  easy  control  over  the 
mean  flow  and  hence  cannot  test  for  these  effects. 

These  advantages  offer  potentially  significant  improvements  in  the  meas¬ 
urement  of  velocity  response  properties  over  the  T-burner  approach.  The 
studies  described  in  the  following  sections  were  conducted  to  assess  the  util¬ 
ity  of  this  new  approach  to  quantitative  velocity  response  measurements. 
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2.2  ANALYTICAL  STUDIES 


Theoretical  analyses  of  the  dual  rotating  valve  were  conducted  as  part  of 
this  program  for  several  purposes.  First,  the  magnitude  and  phase  behavior  of 
the  oscillating  pressure  must  be  predicted  to  see  if  accurate  measurements  of 
the  velocity  response  can  be  made  over  a  range  of  frequencies.  The  sensitiv¬ 
ity  of  the  pressure  oscillations  to  response  function  changes  must  also  be 
estimated.  Finally,  a  method  for  deriving  the  velocity  response  function  from 
the  measured  pressure  oscillations  is  required. 

This  requires  two  formulations  of  the  basic  ballistic  equations  which 
describe  the  oscillatory  motor  behavior.  In  the  first  formulation,  the 
response  of  the  burning  rate  to  velocity  oscillations  and  pressure  oscillations 
is  assumed  to  be  known.  The  equations  are  formulated  to  predict  the  pressure 
and  velocity  oscillations  based  on  these  known  responses.  The  solution  of 
these  equations  provides  a  method  for  estimating  the  behavior  of  the  burner  and 
the  sensitivity  of  this  behavior  to  propellant  changes.  The  second  formulation 
assumes  the  pressure  oscillations  are  known  and  the  velocity  and  pressure 
responses  of  the  propellant  are  unknown.  This  second  formulation,  which  is  the 
inverse  of  the  first  formulation,  is  required  to  derive  the  propellant 
responses  from  experimental  data. 

The  analyses  are  complicated  by  the  fact  that  two  qualitative  models  have 
been  proposed  to  describe  velocity  coupling  a  iinear  model  and  a  non¬ 

linear  model.  In  the  linear  model,  the  oscillatory  burning  rate  responds  to 
the  local  magnitude  and  direction  (relative  to  the  mean  speed)  of  the  velocity 
oscillations.  In  other  words,  the  burning  rate  response  depends  on  the  direc¬ 
tion  of  the  mean  gas  flow  as  well  as  the  magnitude  and  direction  of  the  oscil¬ 
lating  velocity,  but  does  not  depend  on  the  magnitude  of  the  mean  speed.  This 
model  has  been  used  extensively  to  interpret  velocity  coupled  T-burner  data,®’® 
motor  data, ^  and  in  predictions  of  motor  stability.® 

The  nonlinear  model  evolved  by  analogy  to  erosive  burning.  Erosive  burn¬ 
ing  is  postulated  to  result  from  increasing  the  heat  transfer  to  the  burning 
propellant  surface  produced  by  increasing  the  gas  flow  parallel  to  the  surface. 

Q  IQ 

Analytical  studies,  as  well  as  experimental  data,  suggest  there  is  a 
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minimum,  or  threshold,  speed  which  must  be  exceeded  oefore  the  effect  of 
parallel  gas  flow  on  the  burning  rate  becomes  significant. 

Extrapolating  a  threshold  speed  to  oscillatory  flows  is  a  central  feature 
of  the  nonlinear  model.  The  basic  hypothesis  assumes  the  burning  rate  responds 
to  velocity  oscillations  when  the  instantaneous  total  velocity  (mean  velocity 
plus  the  oscillating  velocity)  exceeds  the  threshold  speed.  Analytical  formula¬ 
tions  which  incorporate  this  concept  were  included  in  the  Standard  Stability 
Prediction  Program11  but  the  concepts  have  never  been  tested  experimentally. 

A  critical  review  and  evaluation  of  these  two  models  was  the  subject  of 
a  recent  JANNAF  workshop.12  The  basic  conclusion  from  this  workshop  was  that 
both  models  have  certain  aspects  which  are  reasonable,  but  other  aspects  which 
are  not  acceptable.  Hence,  both  models  were  considered  in  the  analytical 
studies  conducted  as  part  of  this  program. 


2.2.1  Burner  Analysis  Incorporating  Linear  Velocity  Coupling 

2.2. 1.1  Solutions  for  Pressure  and  Velocity  Oscillations  for  Given  Propellant 
Properties 

The  initial  analytical  studies  were  conducted  using  the  one  dimensional 
equations  of  motion  for  the  gases  in  the  combustion  chamber,  a  linear  velocity 
coupling  model,  and  the  ideal  gas  law.  After  linearizing  the  momentum  and 
energy  equations  and  retaining  the  oscillatory  parts,  one  obtains: 


de  3M'  3(M  •  M') 

8Z  dr  3Z 


+  o' 


for  the  momentum  equation  and 


F' 
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(1) 


3f ' 
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+  KH'+ 


3M' 
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,  M*M ' 

+  tv  nr 
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for  the  energy  equation.  The  factor  fi/IMI  arises  in  the  velocity  coupling  term 
because  the  mass  flow  and  flame  temperature  oscillations  are  assumed  to  depend 
on  the  oscillation  in  gas  speed  parallel  to  the  propellant.  A  positive  per¬ 
turbation  in  velocity  gives  a  positive  perturbation  in  speed  if  the  mean  veloc¬ 
ity  is  also  positive.  If  the  mean  velocity  is  negative,  a  positive  velocity 
perturbation  produces  a  negative  perturbation  in  speed.  It  should  be  noted 
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that  these  equations  are  identical  to  those  used  by  Culick1^  to  estimate 
the  linear  stability  of  acoustic  chambers  for  longitudinal  waves.  The 
solutions,  however,  proceed  along  different  paths  because  the  objectives  dif¬ 
fer.  Culick  was  interested  in  calculating  the  change  in  the  complex  eigen¬ 
value  for  the  chamber  to  determine  the  stability  of  self-excited  acoustic 
waves;  thus,  his  solution  uses  approximations  which  are  valid  only  near  reso¬ 
nant  conditions.  In  this  project,  however,  the  primary  interest  is  the  burner 
response  to  driven  flow  oscillations  at  nonresonant  frequencies;  hence  the 
approximations  which  are  appropriate  for  resonant  conditions  are  not  valid  in 
this  development. 


To  obtain  the  solution  for  nonresonant  conditions,  equations  2  and  3  can 
be  combined  to  eliminate  dM'/dZ  by  assuming  separation  of  variables  and 
replacing  3 (  '  J/9 T  by  i  A(  ).  The  result  is  then  differentiated  with  respect  to 
Z.  Terms  containing  M'  and  dM'/dZ  can  then  be  eliminated  to  yield 


dz2 

where 


(v  «  +  )+  v(a  -  q" /d-)]||  -  Tk(0  +  dM/dZ)  l-o  (3) 


Q  =  iX 


1  +  C  +  iXT 
m _ d 

1  +  i\T, 


+  -  rft>  f 


(*») 


Tk=  iX  -  A^qL/S 


(5) 


Y  = 


+  Cl) 
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) 


dM 
*  dZ 


(6) 


Two  approximations  have  been  made  in  the  derivation  of  equation  3.  First, 

M which  limits  the  analysis  to  low  mean  flow  speeds.  Second,  MY<<Q  has 
been  assumed.  If  the  magnitude  of  the  velocity  response  is  of  the  order  of 
unity,  then  fiY  is  of  order  M2.  Q,  on  the  other  hand,  is  of  order  M.  At  high 
magnitudes  of  velocity  response,  this  approximation  breaks  down.  This  point 
will  be  discussed  further  in  section  2.2.2. 3* 


2-6 


i 


Equation  3  is  second  order  and  has  variable  coefficients.  By  letting 


r  =  (CZ  -  N)  /C 


(7) 


c*°-5(i)(Q  +  vf) 

,.^.5T(<Lt^4  |  (o*V’g) 


k  =  Tk(Q  +  dM/dZ) / (4C) 

equation  3  can  be  transformed  to  a  Kummers  equation. 

j2a 


in 


3 -  (MS-* 


The  solution  becomes 


15 


=  D  K  (k,j,r)  +  Er2  K  jk  + 


(8) 

(9) 

(10) 

(ID 


(12) 


The  constants  of  integration  can  be  evaluated  from  the  nozzle  flow  equa¬ 
tions  at  each  end  of  the  chamber.  Assuming  the  nozzles  and  valves  behave  in 
a  quasi-steady  manner,  one  finds 


at  Z  =  0  and 


at  Z  =  1 . 


A  =  M  U  +0.5  (T/T)  \  =  M  T.  £  -  (de/dZ) 

O  o(rO  oj  OKO  O 

Mi  =  Mj  +0.5  (T/T)  =  Mx  T^  - 


(de/dZ) 


(13) 


(in) 


At  this  point,  there  is  a  dilemma  in  relating  the  nozzle  equations  for 
the  oscillating  Mach  number  at  each  end  of  the  burner  to  the  values  derived 
from  equations  1  and  2.  The  nozzle  equations  are  in  terms  of  velocity  and  tem¬ 
perature  oscillations  whereas  equations  1  and  2  are  written  in  terms  of  velocity 
and  pressure  oscillations.  If  the  flame  temperature  oscillation  were  known 
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separately  from  the  burning  rate  oscillation,  it  would  be  possible  to  introduce 
an  equation  for  the  entropy  oscillation  in  the  burner.  The  solution  to  this 
entropy  equation  would  yield  the  gas  temperature  oscillation  at  the  nozzle 
entrance  in  teras  of  the  pressure  and  the  flame  temperature  oscillations. 
Without  explicit  knowledge  about  the  flame  temperature  oscillation,  however,  an 
approximation  must  be  made.  Conventionally,  the  entropy  oscillation  at  the 
nozzle  entrance  is  assumed  to  be  zero,  which  yields 


With  this  approximation,  the  nozzle  boundary  conditions  become 


at  Z  =  0  and 


at  Z  =  1. 


(^f)  o  +  Vl  •  Q"o  [°'5fo  <7-  »  +*o] 

(^i)  +V/i  ■  <*1  [°-5'i  <7-  +*i] 

1 


(15) 


(16) 


(17) 


Because  Y  changes  sign  when  the  mean  flow  Mach  number  goes  to  zero,  two 
sets  of  constants  are  required:  one  set  for  0<Z<Xi  and  one  forX1<Z<1.  The 
additional  boundary  conditions  are  obtained  by  requiring  the  oscillatory  pres¬ 
sure  and  velocity  to  be  continuous  for  Z  s  X.,. 

Thus,  if  the  acoustic  response,  particle  damping  parameters  and  chamber 
dimensions  are  known,  the  oscillating  components  of  pressure  and  velocity  can 
be  predicted  for  any  point  in  the  chamber.  Furthermore,  these  predictions  can 
be  made  for  any  phase  relationship  between  the  two  values,  simply  by  selecting 

A  A 

appropriate  values  of  4/q  and  y/j. 

There  are  several  significant  features  to  this  solution.  First,  momentum 
and  energy  effects  are  included  along  with  the  transient  mass  balance  contribu¬ 
tions.  Second,  the  analysis  includes  pressure  coupled  response,  flow  turning 
and  particle  damping  effects.  Third,  the  analysis  permits  the  study  of  the 
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two-valve  configurations  as  well  as  the  one-valve  configuration.  Fourth,  the 
analysis  is  not  restricted  to  low  frequencies.  Thus,  equation  12  is  a  general 
solution  which  should  be  valid  for  one-dimensional  burners  at  all  frequencies. 
This  model  was  programmed  in  Fortran  for  the  Burroughs  6700  computer.  The 
basic  definition  for  the  Kummers  function1^  was  used  to  calculate  these 
functions. 

2.2. 1.2  Verification  of  Solution 

The  solution  of  the  equations  described  in  the  preceding  paragraphs  was 
verified  in  two  ways.  First,  at  low  frequencies,  the  equations  should  reduce 
analytically  to  the  equations  derived  for  the  pressure  coupled  rotating 
valve,  » 17  when  the  velocity  response  is  set  at  zero.  For  low  frequencies, 
the  Kummers  functions  in  equation  12  approach  unity.  With  this  simplifica¬ 
tion,  equation  12  leads  to  the  identical  result  reported  in  references  16  and 
17  for  the  case  of  one  valve  and  for  both  valves  when  operated  in  phase. 
Numerical  calculations  also  yield  identical  results.  This  is  an  important 
point  since  it  provides  a  numerical  verification  of  the  computer  programming 
in  addition  to  the  analytical  verification. 

Second,  this  model  should  predict  the  correct  acoustic  damping  when  the 
combustor  is  driven  at  frequencies  near  the  natural  acoustic  frequencies.  In 
particular,  one  would  expect  the  frequency  difference  at  the  half-power  ampli¬ 
tudes  (i.e.,  0.707  x  the  peak  amplitudes)  would  be  related  to  the  overall  sys¬ 
tem  damping  of  the  self-excited  system  by  the  expression 

a  jAf  (18) 

f  '  f 

The  left  side  can  be  evaluated  independently  from  Culick’s  solutions 
while  the  right  side  can  be  evaluated  from  numerical  solutions  of  the  model. 

This  comparison  has  been  made  in  three  cases  and  the  results  are  shown  in 
table  2-1 . 

The  first  case  contained  only  pressure  coupling  effects  and  used  one 
valve.  The  second  and  third  cases  incorporated  both  pressure  and  velocity 
coupling,  as  well  as  particle  damping  effects.  In  case  two,  the  response 


TABLE  2-1.  COMPUTED  SYSTEM  DAMPING  (<*/f) 


TM 126 


Case _ Culick  Analysis _ Equation  (12) 

1  0.112  0.117 

2  0.155  0.156 

3  0.102  0.108 


functions  were  low  (i.e.,  0.2)  while  in  case  three,  they  were  approximately  an 
order  of  magnitude  higher. 

The  excellent  agreement  found  between  these  two  methods  substantiates  the 
analysis.  This  demonstrates  the  analysis  is  not  limited  to  the  low  frequency 
bulk  flow  conditions  of  the  rotating  valve,  but  is  also  applicable  to  the 
higher  frequency  conditions  near  the  acoustic  mode  frequency.  Thus,  this 
model  could  provide  a  basis  for  interpreting  data  obtained  in  the  modulated 
throat  motor1®  as  well. 

2.2. 1.3  Solutions  for  Velocity  Response  from  Measured  Pressure  Oscillations 

A  /\ 

The  analysis  in  the  preceding  section  calculates  eQ  and  t 1 ,  for  input 
values  of  the  propellant  responses.  In  combustion  tests  however,  £0  and  ^ 
will  be  measured  and  a  procedure  is  needed  to  derive  the  responses  from  the 
data.  Efforts  to  obtain  analytical  solutions  were  not  successful.  Since  the 
responses  are  implicit  in  equation  12,  rearranging  the  forward  solution  to 
obtain  expressions  which  are  explicit  in  the  response  was  also  unsuccessful. 

This  left  two  alternatives:  a  trial  and  error  approach,  and  an  approxi¬ 
mate  solution.  In  the  trial  and  error  approach,  one  would  assume  the  response 
value  and  estimate  the  expected  magnitude  and  phase  for  the  pressure  oscil¬ 
lations.  These  would  be  compared  to  the  data;  new  values  for  the  response 
would  be  estimated;  and  new  values  for  the  pressure  oscillations  would  be  cal¬ 
culated.  Attempts  to  apply  this  approach  met  with  little  success.  The  prin¬ 
cipal  difficulty  was  in  developing  a  method  for  estimating  the  new  response 
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value  from  the  comparison  of  the  estimated  and  the  observed  pressure  oscilla¬ 
tions  . 


Attention  was  then  directed  to  developing  an  approximate  solution.  A 
large  number  of  numerical  calculations  were  made  to  examine  the  predicted 
behavior  of  the  pressure  and  velocity  oscillations.  To  conduct  these 
calculations,  a  method  of  selecting  response  functions  was  needed.  A  rational 
basis  is  provided  by  using  thermal  wave  combustion  models. 19»20»21 


r  =  n(l  +  A  -  R) 

p  +  A/\l  -  R 


(19) 


and 


r  =  m(1  +  A  -  R) 
v  +  A/X1  -  R 


(20) 


The  combustion  parameters,  A  and  R,  can  then  be  determined  from  pressure 
response  measurements.  The  parameter,  m,  must  be  selected  arbitrarily,  although 
unity  appears  to  give  reasonable  values,  based  on  limited  experimental  values 
of  velocity  response  functions  from  T-burner  tests. 


Using  values  for  the  aluminized  propellant,  ANB  3066,  (A  =  15,  R  =  ^ . 97 , 
n  =  0.27,  m  =  1 ,  T=  31.5)  and  the  valves  out  of  phase  by  180°,  a  large  number 
of  calculations  were  run  for  the  experimental  geometry  to  be  used  in  this 
study.  These  calculations  showed  that  at  low  frequencies  (i.e.,  A=  0.75)  the 
oscillating  velocity  is  essentially  spatially  uniform.  The  magnitude  is  deter¬ 
mined  strictly  by  the  size  of  the  oscillating  area  of  the  two  valves  relative 
to  the  time-averaged  area.  In  addition,  the  oscillating  pressure  is  shown  to 
be  nearly  a  linear  function  of  axial  position. 


Using  the  linear  pressure  observation,  an  approximate  solution  of  the 
energy  equation,  which  is  explicit  in  the  velocity  response,  can  be  developed. 


(R  +  co  ) 

V  V 
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The  accuracy  of  this  approximation  has  been  examined  numerically.  The 
exact  solution  was  used  to  predict  oscillating  pressures  and  velocities  using 
estimated  combustion  parameters  for  ANB  3066.  These  pressures  and  velocities 
were  then  used  to  simulate  experimental  observations  as  inputs  to  equation  21. 
The  velocity  response  derived  from  equation  21  was  then  compared  with  the  orig¬ 
inal  input  value  to  determine  the  accuracy  of  the  approximate  solution.  Fig¬ 
ures  2-3  and  2-h  show  comparisons  of  the  derived  and  input  velocity  responses 
for  a  variety  of  conditions.  This  comparison  was  made  using  the  imaginary 
part  of  the  response  since  this  is  a  parameter  used  in  motor  stability 
predictions. 

The  results  in  figure  2-3  were  obtained  by  assuming  the  propellant  con¬ 
tained  no  metal.  With  the  valves  phased  by  exactly  180°,  the  maximum  error 
between  the  input  and  derived  response  is  ±10$  at  the  point  where  the  response 
is  a  maximum.  Allowing  for  phase  misalignments  of  3*5°  between  the  two  rota¬ 
ting  valves,  this  error  increases  to  ±20$.  This  error  is  typical  for  the 
experimental  conditions  and  hence  ±20$  accuracy  is  expected  for  this  procedure. 

Figure  2-H  shows  the  effect  of  particle  damping  on  the  comparison  between 
the  exact  and  derived  velocity  response.  These  calculations  were  made  by 
using  the  maximum  particle  damping  at  each  frequency.  Again,  the  agreement  is 
better  than  ±10$.  The  approximate  analysis,  therefore,  appears  to  provide  a 
reasonably  accurate  solution  at  the  low?.'  frequencies  for  the  conditions 
examined.  This  approximate  solution  also  suggests  a  good  method  for  estimat¬ 
ing  initial  inputs  if  the  complete  trial  and  error  analysis  is  ever  required 
to  improve  the  accuracy  of  the  data  reduction  process. 

2. 2.1.1*  Second  Order  Acoustic  Analysis 

Experimental  studies,  which  are  described  in  detail  in  section  2.3» 
showed  significant  generation  of  frequencies  which  are  harmonics  of  the  driven 
frequency.  Since  the  waveform  of  the  area  oscillations  has  low  harmonic  con¬ 
tent,  ^  they  must  be  generated  by  nonlinear  burner  acoustics  or  by  nonlinear 
combustion  processes.  Therefore,  an  analysis  was  conducted  to  estimate  the 
harmonic  content  produced  by  second  order  acoustics  and  linear  velocity 
coupling. 
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Figure  2-*J.  Error  Analysis  Study,  ANB  3066  (Particle  Effects  Included) 
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In  this  second  order  analysis,  each  variable  is  written  in  the  form: 

e=e°  +  e'  +  e" 

By  expressing  velocity,  density,  burning  rate,  and  flame  temperature  os¬ 
cillations  in  this  manner,  separate  conservation  equations  can  be  derived  for 
the  zeroth,  first  and  second  order  perturbations.  The  second  order  equations 


.  3M" 

dr 


3(M  .  M") 


+  —  .  M"  =  -  — 
dZ  2 


1  3(M’)‘ 

2  dZ 


3(M  H') 


dM°  m'bM’ 


de  "  3M"  dM°  " 

dr  dZ  7  az  e 


+  M 


°  djT 

dz 


II 

o|T 

dM  _f_ 
dZ  T  0 


+ 


(2H) 


del  ,  dM° 
3Z  +  dZ 


In  these  equations,  terms  of  order  (M°) 2  have  been  neglected  and  the 
density  perturbation  is  related  to  the  pressure  oscillation  using  the  isen- 
tropic  approximation  . 


There  is  a  problem  defining  ^  and  ^-1  .  In  the  first  order  equations, 
only  the  sum  of  these  terms  occurs,  as  it  does  on  the  left  hand  side  of  equa¬ 
tion  24.  These  terms,  however,  also  appear  separately  in  equation  24.  To  sep¬ 
arate  ^  from  ^  it  was  assumed  that  ^  is  related  isentropically  to  e ' . 
Corrected  first  order  response  function  values  were  then  used  for  . 


The  solution  for  the  second  order  perturbations  was  obtained  by 
expressing  the  equations  in  finite  difference  form.  The  burner  was  divided 
into  10  stations  and  the  equations  were  written  for  the  segments  between  the 
stations.  This  gave  20  equations  for  the  variables  c" i  through  f"io  and  M"i 
through  M" iq  (nine  momentum  equations,  nine  energy  equations  and  a  boundary 
condition  at  each  end).  The  first  order  solution  is  also  needed  at  each  sta¬ 
tion.  Since  the  left  hand  sides  of  the  first  and  second  order  equations  are 
identical,  the  first  order  solution  was  found  in  the  same  manner.  These 
values  were  then  used  in  the  right  hand  side  of  the  second  order  solution. 

The  product  of  two  first  order  perturbations  gives  a  steady  component  and  a 
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component  at  double  the  frequency  of  the  first  order  perturbation.  Friedly^2 
showed  that  if  e*  =  Re|c  e^^J  then 

(  e')2  *  0.5  Re  ^  (  e  )  (  e  )*  +  (e  )2e12Wt 

The  results  of  this  analysis  for  the  same  cases  treated  in  the  parametric 
study  in  section  2.2. 2. 2  are  shown  in  table  2-2.  The  double  frequency  ampli¬ 
tudes  produced  by  second  order  acoustics  are  clearly  very  small.  Hence, 
nonlinear  acoustics  are  not  expected  to  be  a  significant  factor  in  the 
generation  of  harmonic  frequency  components. 

2.2.2  Burner  Analysis  Incorporating  Nonlinear  Velocity  Coupling 
2.2.2. 1  Formulation  of  Nonlinear  Velocity  Coupling  Model 

The  nonlinear  velocity  coupling  model  used  in  this  analysis  combines  some 
kinematic  concepts  extrapolated  from  steady  state  erosive  burning  with  a  fre¬ 
quency  dependent  response  function.  This  model,  which  has  been  studied  by 
Price,2, ^  Culick14  and  Condon, 23  was  developed  by  heuristic  arguments  to 
qualitatively  explain  some  observed  instabilities. 


The  basic  approach  assumes  the  velocity  oscillations  generate  an  oscilla¬ 
tory  heat  flux  to  the  burning  propellant  surface.  Since  the  heat  flux  depends 
on  the  magnitude  of  the  velocity  and  not  the  direction,  the  oscillatory  heat 


TABLE  2-2.  PRESSURE  AMPLITUDES  FOR  LINEAR  COUPLING  WITH  NONLINEAR  ACOUSTICS, 

ANB  3066,  x  =  0.1256 

1  T6194 

Frequency _ At  Driver  Frequency  At  Double  Frequency _ Ratio,  dB 


0.0598 

0.616*2 

0.998-5 

-56 

0.1197 

0.484-2 

0.765-5 

-56 

0.1795 

0.438"2 

0.636-5 

-57 

0.239*1 

0.412"2 

0.562-5 

-57 

0.2992 

0.37 2~2 

0.560-5 

-56 

0.3590 

0.298-2 

0.567*5 

-54 

0.4189 

0.206-2 

0.574-5 

-51 
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flux  to  the  surface  is  rectified.  A  threshold  speed  is  then  added  to  the 
model  by  requiring  the  instantaneous  speed  to  exceed  the  threshold  before 
oscillatory  effects  are  produced  in  the  burning  rate.  These  basic  concepts  are 
shown  in  figure  2-5.  The  driving  velocity  oscillation  is  compared  in  fig¬ 
ure  2-5A  to  the  mean  velocity  and  the  threshold  speed.  The  "effective"  velocity 
(i.e.,  the  heat  flux)  is  shown  in  figure  2-5B.  Mathematically,  tuis  model  can 
be  formulated  as  follows.  The  effective  oscillatory  velocity  can  be  written: 


Figure  2-5.  Effective  Velocity  for  Nonlinear  Velocity  Coupling 
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In  equation  25,  h(  )  is  the  Heaviside  unit  step  function,  which  is  unity  if  the 
argument  is  positive  and  zero  if  the  argument  is  negative. 


The  term  is  the  amplitude  of  the  rectified  wave  at  the  driver  fre¬ 
quency.  This  term  is  then  combined  with  a  "response  function"  to  produce 
burning  rate  (and  flame  temperature)  oscillations  according  to  the  equation 


A  A 

rab  .  Tf 


Z^+Z1  =  <Rd+<V  7e  +  <Rv+Wv)*l 


A  more  complete  analysis  of  this  model  can  be  found  in  papers  by  Price 
and  Dehority,^'^  Culick,^  and  Condon, ^  and  in  the  Standard  Stability  Predic¬ 
tion  Program  Manual.11 

A  _ 

For  given  values  of  M,  M  and  Mth  relatively  simple  expressions  can _ 

.  .  _  |M|  |M| 

be  derived  for  /^/M  which  are  valid  for  a  particular  domain  of  —  and  . 

Figure  2-6  shows  a  plot  of  iftl/M^h  versus  IMf/M^  in  which  four  distinct 
regions  are  identified: 


In  region  A: 


Zi  =  0  ( iMl  +  IMI  <  Mth) 

M 


In  region  B: 


In  region  C: 


tlmJL  (  |m!  -  IMI  >  M  ) 
m  I'm  l 


-  m  , )  cos  e ,  +  -M"  v  2ei  +  sin  <2e 

tn  i  1 


1  J  IMl 


for  ( I M I  -  |M|  <  M  h) 


and  (  I M I  +  |M|  >  M^) 
and  (  |M|  -  |M|  <  Mth> 
where  6^  =  sin  1  j(Mth  “  I M |  )  /  |M|  J 
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In  region  D: 


tl  -  2  [(  |M!  -  Mth)  cos 

A  >r|A|  l 

♦  2  [(  |M|  +  Mth)  cos  02 

fflAi  [ 


6\  ♦  |  Mi  (  v  -  2  9]  *  sin  ( 2  &  i  )~|  -21- 

14  -I  |M| 

_  I M  (  „  -  202  +  sin  (20?)]_M 

11  J  |RI 


for  (  |M|  -  |M|  >  M  ) 
tn 

0j  =  sin  * 

[(Mth  -  |M|  )/|M|J 

a  •  -1 

a  2  sin 

(Mfch  +  IM|  )  /  (M|l 

(3D 


q  i  i  i 
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Figure  2-6.  Domain  of  Various  Types  of  Velocity  Coupling 

21058 


2-18 


These  formulas  were  coded  for  a  digital  computer.  The  validity  of  the 
calculations  was  established  by  comparing  these  results  with  direct  computa¬ 
tions  of  the  Fourier  coefficients  from  the  time  domain  wave  form. 

A  similar  analysis  was  performed  to  obtain  and  /jg*  The  ^2  is 
interest  because  it  produces  a  harmonic  of  the  driven  frequency  and  the  burner 
data  show  significant  harmonic  content. 


In  regions  A  and  B: 


In  region  C: 


“o  ■  0 
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In  region  D: 
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In  regions  A  and  B: 


Plots  of  p.,/M  and  p are  shown  in  figures  2-7  and  2-8.  The  figures 
shew  that  the  fundamental  component  of  the  rectified  velocity  varies  almost 
linearly  with  M  for  a  constant  amplitude  of  the  driving  velocity  oscillation. 
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The  second  component  depends  on  both  M  and  M  in  a  manner  that  is  not  easily 
generalized. 

2. 2. 2. 2  Solution  of  Ballistic  Equations  Incorporating  Nonlinear  Velocity 
Coupling 

The  nonlinear  velocity  coupling  model  can  now  be  incorporated  into  the 
ballistic  equations  for  the  burner.  Essentially  this  amounts  to  replacing  the 
linear  velocity  coupling  terms  in  equation  2  with  nonlinear  velocity  coupling 
terms.  This  presents  an  additional  complication.  The  nonlinear  terms  vary 
with  axial  position  which  changes  the  form  of  the  second  order  wave  equation 
derived  fran  equations  1  and  2. 


This  effect  can  best  be  explained  by  examining  the  velocity  distributions 
in  the  burner  as  shown  in  figure  2-9.  This  plot  shows  M  versus  Z.  The  locus 
of  the  minimum  and  maximum  values  of  M(Z)  during  oscillation  are  shown  dotted 
above  and  below  M.  Two  horizontal  lines  representing  the  plus  and  minus  thres¬ 
holds  are  also  shown.  The  burner  is  divided  into  five  regions  since  each 


Mean  plus 
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Figure  2-9.  Description  of  Regions  Containing  Nonlinear  Velocity  Coupling 
V-10H69  21061 
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region  has  a  different  formulation  for  the  velocity  coupled  driving.  In 
regions  1  and  5,  linear  velocity  occurs  since  the  M  +  M'  is  always  greater  than 
the  threshold  value.  In  regions  2  and  H  nonlinear  velocity  coupling  occurs, 
while  no  response  is  generated  in  region  3- 

This  situation  requires  the  solutions  of  equations  1  and  2  for  each 
region  of  the  burner.  The  linear  solutions  developed  in  section  2. 2. 1.1  can  be 
used  directly  in  regions  1  and  5,  and  in  region  3  by  setting  the  velocity 
response  to  zero.  Regions  2  and  H  are  more  difficult  since  both  M  and  M  vary 
with  axial  position. 


One  could  always  resort  tc  numerical  integrations  which  incorporate  these 
effects  at  each  axial  position  in  regions  2  and  H.  However,  the  basic  objec¬ 
tive  of  this  analysis  is  to  predict  the  low  frequency  behavior;  figure  2-7  sug- 

S\ 

gests  an  approximate  approach.  At  low  frequencies,  M  is  nearly  independent  of 
axial  position.  Since  M  is  linear  with  axial  position,  figure  2-7  suggests 

/.  /s 

/u  /m  is  a  linear  function  of  axial  position  in  regions  2  and  4. 

By  using  this  approximation,  a  simple  expression  for  the  velocity  coupled 
driving  in  each  of  the  five  regions  can  be  written  by  redefining  the  velocity 
coupled  driving  term,  Y.  For  nonlinear  velocity  coupling,  Y  becomes 


dM 

Y  *  -  77  (R  +  w) 

dZ  v  v 


in  region  1 


(38) 


Y  = 


£«,♦“*>  Ki  r-rr. 


'1 
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(39) 


Y  =  0 


in  region  3 


(HO) 


dM  ^  “  ^3 

Y  =  ^  (Rv  +wv)  y—  z  in  region  4 

4  3 


(HI) 


dM 

Y  -  5?  (R  +  w  ) 
dZ  v  v 


in  region  5 


(U2) 
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The  second  order  differential  equation  11  can  be  derived  if  the  constants 
C  and  N  are  redefined  as 


_  ,  _  -  dM  ,  1  dM  „ 

T  +  Q  +  2—  +  -—  (R+oj)  G 
k  dZ  Q  dZ  v  v  c 


(R  +w  )  ^  G 
v  v  dZ  c 


dM 

dZ 


(43) 


2\  =  X  — 
ldZ 


T*+Q  +  2§  +  ^§  (Rv+cV 


(R  +W  )^G 
v  v  dZ  n 


dM 

dZ 


(44) 


where 


G  =  -K  /(Z  -  Z  )  and  G  =  Z  G  in  region  2 

c  11  2  n2c 

G  =  K./(Z.  -  Z.)  and  G  =  Z.G  in  region  4 
c  4  4  3  n  3  c  6 


at  Z , 


(45) 

(46) 


(47) 


K 


4 


at  Z . 


(48) 


The  constants,  K-j  and  ,  ar  1  if  regions  1  and  5  exist.  There  are 
conditions,  however,  where  these  regions  do  not  exist  and  hence  where  or 
will  be  between  0  and  1. 

The  solution  is  essentially  the  same  as  in  section  2. 2. 1.1,  except 
there  are  five  regions  instead  of  two  as  in  the  linear  case.  Some  iteration 
is  required,  however,  since  the  locations  of  the  boundaries  between  regions 
depend  on  M(Z),  which  is  rot  known  until  the  final  solution  is  known. 
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As  a  first  guess,  M  at  each  end  of  the  burner  is  calculated  from  equations 
16  and  17  with  e  =  0.  Intermediate  values  are  assumed  to  be  linear  between 
the  ends.  From  these  estimates  of  M,  the  locations  of  the  region  boundaries 
can  be  estimated.  The  solution  in  each  region  is  defined  in  terms  of  the 
Rummers  functions.  Equations  are  written  equating  e  and  M  at  the  interfaces 

*  A 

between  the  regions  together  with  the  boundary  conditions  for  t  and  M  at 
Z  =  0  and  Z  =  1.  A  total  of  10  linear  complex  equations  and  10  unknowns 
can  be  solved  to  obtain  the  coefficients  D  and  E  in  each  region.  From 
these  coefficients,  M  at  the  interfaces  can  be  established  and  new  estimates 
for  the  locations  of  the  interfaces  can  be  made.  This  iteration  is  continued 
until  the  changes  in  the  locations  of  the  interfaces  are  as  small  as  desired. 

In  practice  only,  one  iteration  is  usually  required. 

A  FORTRAN  computer  code  was  written  for  this  analysis  on  the  Burroughs 
6700.  Besides  the  normal  checking  process,  cases  were  run  with  a  zero 
threshold  speed  and  low  amplitude  area  oscillations.  Under  these  conditions 
the  nonlinear  model  becomes  identical  to  the  linear  model  (except  for  tie 

A  __ 

small  region  where  I  Ml  >  M -  Excellent  agreement  with  the  linear  model 
was  observed. 

A  parametric  study  was  then  performed  assuming  the  propellant  responses 
followed  the  thermal  lag  model.  Figures  2-10  through  2-12  show  the  results 
for  a  typical  propellant  and  burner  geometry.  In  figure  2-10  the  amplitude  of 
pressure  oscillation  follows  the  general  trend  of  decreasing  amplitude  with 
frequency  expected  for  the  bulk  mode.  Figure  2-11  shows  the  ratio  of  the  pres¬ 
sure  amplitudes  at  each  end  of  the  burner.  The  corresponding  phase  angles  are 
shovm  in  figure  2-12.  In  the  amplitude  ratio  vs  frequency  curve,  the  minimum 
tends  to  shift  to  a  lower  frequency  as  the  threshold  velocity  increases.  The 
plot  of  phase  in  figure  2-12  shows  that  the  Mth  =  0.03  curve  has  a  different 
character  from  the  other  three  curves.  For  =  0.03,  all  the  velocities  are 
less  than  the  threshold  and  nc  coupling  occurs.  There  is  a  drastic  change 
frcm  Mth  =  0.02  to  0.03.  Additional  cases  were  run  to  study  this  change  in 
greater  detail.  These  cases,  which  are  presented  in  figure  2-13,  show  that  as 
the  amount  of  velocity  coupling  approaches  zero,  the  character  of  the  phase 
curve  changes  rapidly. 
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Figure  2-10.  Pressure  Amplitude  from  Analytical  Model  of  Velocity-Coupled 
Rotating  Valve  (Input  Response  Based  on  ANB  3066) 
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Figure  2-11.  Ratio  of  Pressure  Amplitudes  (Input  Response  Based  on  ANB  3066) 
V-101169  2106 
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Figure  2-12.  Phase  of  Pressure  at  Z  =  0  to  Pressure  at  Z  =  1  (Input  Responses 

Based  on  ANB  3066)  21064 
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Figure  2-13.  Effect  of  Mth  cn  Phase  of f  vs  Z 
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Approximate  solutions  for  amplitude  of  the  pressure  oscillations 
generated  at  the  higher  frequency  components  were  also  made  using  this  anal¬ 
ysis.  The  oscillations  in  gas  speed  can  be  decomposed  into  the  Fourier  compo¬ 
nents  and  each  component  can  be  assumed  to  drive  the  propellant  independently. 
Once  the  first  order  solution  is  obtained,  the  burn  rate  oscillation  at  the 
double  frequency,  for  example,  can  be  obtained  by  integrating  (Rv  +  v ) DF ^ 2 
over  the  regions  containing  nonlinear  velocity  ooupling.  Assuming  the  burner 
operates  in  the  bulk  mode,  the  pressure  oscillation  at  the  harmonic  frequency 
can  be  estimated  from  the  equation 


-O’v+y  prfc) 

7  (R  +  co  )  +  1-37 

P  P  DF  2 


2iX 

(dM/dZ) 


(149) 


Estimates  for  the  pressure  amplitudes  at  the  second  harmonic  of  the 
driver  frequency  were  made  assuming  combustion  parameters  for  ANB  3066.  These 
results  are  tabulated  in  table  2-3  below. 


TABLE  2-3- 


PREDICTED  PRESSURE  AMPLITUDES  FOR  NONLINEAR  VELOCITY  COUPLING 
WITH  LINEAR  ACOUSTICS  (Mth  =  0.0,  \  =  0.1256) 

1  T6195 


Driver 

Frequency 


Pressure 
Amplitude  of 
Driver 
Frequency 


Pressure 

Amplitude 

Second 

Harmonic 


0.0598 

0.762 

X 

10-2 

0.367  x  10*3 

0.1197 

0.527 

X 

io-2 

0.346  x  10-3 

0.1795 

0.452 

X 

10-2 

0.314  x  10-3 

0.23914 

0.418 

X 

10~2 

0.202  x  10-3 

0.2992 

0.377 

X 

10-2 

0.123  x  10-3 

0.3590 

0.303 

X 

10-2 

0.81  x  10"^ 

0.4189 

0.214 

X 

10-2 

0.59  x  IQ"14 

2- 

-28 

Mr**1 


These  results  show  that  measurable  harmonic  generation  would  be  produced 
by  nonlinear  velocity  coupling  effects.  Since  the  threshold  speed  was  zero  in 
this  case,  the  harmonics  resulted  from  flow  reversal  alone.  This  point  will 
be  discussed  further  in  section  2.3-3  when  harmonic  content  of  observed  pres¬ 
sures  are  compared  to  these  predictions. 


2. 2. 2. 3  Solution  for  Nonlinear  Velocity  Response  from  Measured  Pressure 
Oscillations 

The  approach  to  deriving  the  nonlinear  response  function,  (Rv  +ojv)  from 
measured  pressure  oscillations  is  similar  to  the  approach  used  in  the  non¬ 
linear  case,  section  2. 2. 1.3-  The  energy  equation  for  the  nonlinear  case  can 
be  written 


dM 

dZ 


i7(RP 


+  cop)  -  i  A  -7 


e  -  M  — - 
dZ 


+  cu  ) 

V 


M(Z)  and  e  (Z)  are  assumed  to  vary  linearly  with  axial  position.  Hence,  the 
pressure  can  be  defined  in  terms  of  the  measured  pressure  at  each  end  of  the 
burner  and  the  flov.’  oscillation  can  be  derived  from  the  transient  nozzle  flow 
equations,  i.e.,  equations  13  and  1H.  Equation  50  can  then  be  integrated  to 
yield 
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Since  M(1)  and  M(0)  are  calculated  and  t(0)  and  e(1)  are  measured,  the 
three  integrals  can  be  evaluated  and  a  value  can  be  obtained  for  (Rv  +  wv)est. 

/*A  A  A 

PjdZ  f  a  linear  relationship  between  fjf  and  M  is  assumed 
for  regions  2  and  *4. 


An  effort  was  made  to  improve  the  accuracy  of  this  estimate  by  using  bet- 

A  A  ^ 

ter  estimations  for  e(Z)  and  M(Z).  Equation  51  can  be  used  to  get  a  new  M(Z) 
by  integrating  from  0  to  Z  instead  of  from  0  to  1.  Similarly,  the  momentum 
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equation  can  be  used  to  improve  the  estimate  for  e  (Z). 
written: 


df  -  -  d (M  M) 

dz  =  F  -°-~dz~ 


iXM 


This  equation  can  be 

(52) 


With  an  estimate  for  M(Z)  this  equation  can  be  integrated  to  get  a  new  esti¬ 
mate  of  ciZ).  This  process  could  be  repeated  many  times  to  get  better  esti¬ 
mates  for  (Rv  +  %)est’ 


This  scheme  of  successive  approximation  was  implemented  numerically  by 
dividing  the  chamber  into  50  elements.  The  method  converges  rapidly;  five 
decimal  place  accuracy  is  achieved  within  five  trials. 


The  applicability  of  this  inverse  solution  was  examined  over  a  range  of 
conditions  by  using  the  forward  solution  to  estimate  the  pressure  oscillation 
characteristics .  In  many  cases,  good  agreement  was  found  between  the  derived 
and  input  values  for  the  velocity  coupled  response  function.  In  other  cases, 
substantial  inaccuracies  were  noted,  particularily  at  frequencies  where  the 
magnitude  of  the  velocity  response  is  large.  Figure  2-HJ  shows  examples  of 
both  results.  With  lower  values  of  dM/dZ,  good  agreement  was  found  while 
higher  values  of  dM/dZ  show  significant  discrepancies. 

One  possible  source  for  this  difficulty  is  the  assumption  that  Q»MY 
which  was  carried  over  from  the  linear  velocity  coupling  analysis.  Parametric 
studies  appear  to  show  that  cases  when  this  assumption  is  valid  provide  better 
agreement.  Additional  study  of  this  point  is  needed  to  develop  a  better  under¬ 
standing  of  the  limitations  in  this  data  reduction  procedure. 

So  far,  this  method  has  assumed  the  threshold  speed  is  known.  Consider¬ 
able  effort  was  devoted  to  the  development  of  a  method  for  deriving  this 
threshold  speed  from  the  data.  No  direct  method  could  be  found  for  deriving 
the  threshold  explicitly  from  the  equations.  Parametric  calculations  showed 
that  if  several  values  of  are  used  for  a  given  set  of  t0  and  e1 ,  a  curve 
of  Im(Rv  +  wy)  versus  Mth  can  be  drawn.  Figure  2-15  shows  two  such  curves 
for  different  values  of  x1 .  In  the  forward  solution,  each  set  of  Im(Rv  +  <-cv) 
and  Mth  (for  constant  produces  t(0)  and  e(1)  values  which  are  nearly 
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identical.  Thus,  there  is  really  no  unique  solution  for  and  Im(Rv  +  tov). 
Qualitatively,  this  is  equivalent  to  stating  that  conditions  of  low  response 
and  low  threshold  speeds  cannot  be  separated  from  conditions  of  high  response 
and  high  threshold  speed.  In  other  words,  figure  2-9  shows  a  little  response 
from  a  large  surface  is  equivalent  with  a  large  response  from  a  snaller  sur¬ 
face.  Since  the  burner  is  operating  in  the  bulk  mode,  perhaps  this  result 
should  have  been  expected. 


The  effect  of  errors  in  pressure  measurements  on  the  derived  value  of 
(Rv  +  ^v)  was  also  evaluated  by  using  this  inverse  analysis.  One  of  the  cases 
used  in  the  parametric  study  in  section  2. 2. 2.2  was  used  as  a  base  ease.  Two 
other  cases  were  run:  one  with  e0  changed  by  5 %  and  held  constant,  and 
another  with  changed  by  5%  and  c0  held  constant.  The  error  in  the  response 
produced  by  errors  in  the  pressure  measurements  can  be  written 
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This  equation  can  be  normalized  by  the  base  case  parameters. 
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The  calculated  results  for  X=  0.359  and  Mt^  =  0.01  are: 

.  0.52.’1  0.365e13-005p4L| 


(R 


This  indicates  that  the  percent  of  error  in  (Ry  wy)  can  be  expected  to 
be  about  the  sum  of  the  percentage  error  in  e0  and  .  Similar  results  were 
found  for  other  cases  evaluated.  Hence,  the  data  reduction  method  does  not 
magnify  the  source  of  error. 
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Finally,  the  effect  of  errors  in  the  phasing  between  the  area  oscilla¬ 
tions  at  the  two  ends  was  investigated.  Figure  2-16  shows  the  effect  of  changes 
in  the  phase  angle  on  the  amplitude  of  the  pressure  oscillation,  |£.,|. 

When  combined  with  equation  5*1  this  calculation  shows  the  importance  of 
minimizing  errors  in  phase  angle  to  achieving  reasonable  accuracy  in  the 
derived  response. 

Significant  effort  was  also  devoted  an  alternative  approach  for  deriving 
the  velocity  response  from  experimental  data.  This  approach  is  based  on 
Newton's  Method.  The  Rummers  function  analysis  in  section  2. 2. 1.1  calculates 
eo  ^p  ♦  wpi  Rv  +  cjv*  and  (Rp  +  cjp,  Ry  +  wy,  Mt^).  The  effect 

of  small  changes  in  (Rp  +  cop)  and  (Ry  +  on  eD  and  can  be  expressed 
as: 


0  0  1  0.2  0  3  0  4  0  5  0  6 
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Figure  2-16.  Effect  of  Phasing  of  Driving  on  Amplitude  of  6i 
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These  equations  can  be  used  to  set  up  an  iterative  procedure  for  obtaining 
(Rp  +  Wp)  and  (Rv  ♦  wv)  given  e0 ,  ei  and  Mth 
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(R  +  co  ).  is  the  i  th  estimate  of  (R  +  w  ) 

P  P  1  P  P 

(R  +  u>  )  .is  the  i  th  estimate  of  (R  +co  ) 
v  v  i  v  v 

f  .  is  the  value  of  £  obtained  in  the  Rummers  function  analysis  using 
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.  is  the  value  of  obtained  in  the  Rummers  function  analysis  using 
(R  +oJ  ).  and  (R  +  oJ  ). 

p  pi  V  VI 

i  is  the  measured  value  of  i 

Co  cO 

is  the  measured  value  of 


These  equations  can  be  used  to  solve  for  (Rp  ♦  Wp)^  and  (Rv  + 

The  partial  derivatives  must  be  calculated  from  the  analysis  by  incrementing 
the  independent  variables  slightly. 
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This  procedure  worked  using  the  e0  and  from  the  forward  analysis  to 

recover  the  values  of  (Rp  +  wp)  and  (Rv  +  wy) ,  when  the  value  of  was  con¬ 
sidered  known.  However,  when  studies  were  made  to  extend  the  method  to  derive 
values  of  Mth,  the  method  did  not  converge.  Again  the  basic  problem  seems  to 
be  the  inability  to  separate  combustion  conditions  of  low  response  and  low 
threshold  from  conditions  of  high  response  and  high  threshold. 

A  variation  of  this  method  was  also  studied.  The  amplitude  and  phase  of 
the  pressure  oscillation  are  measured  at  each  end  of  the  burner.  With  four 
equations  and  assuming  the  real  part  of  (Rp  +  wp)  is  known  from  a  single  rota¬ 
ting  valve  test,  it  should  be  possible  to  solve  for  Im(Rp  +  wp),  0*v  +tV  and 
This  approach  was  explored  but  again  solutions  could  not  be  obtained. 
Further  study  showed  the  equations  to  be  ill-conditioned  and  development  of 
this  approach  was  terminated. 

2.3  EXPERIMENTAL  STUDIES 
2.3.1  Apparatus  Description 

A  dual  rotating  valve  apparatus  was  designed  and  constructed  as  an  inte¬ 
gral  part  of  this  program.  The  basic  apparatus  layout  is  shown  schematically 
in  figure  2-17.  Figure  2-1 8  shows  the  essential  components  of  the  rotating 
valve  in  exploded  view.  The  rotor  with  a  graphite  sleeve  is  shown  together 
with  four  slots  in  their  relative  positions.  The  two  slots  on  the  upper  left 
of  the  picture  vent  the  two  auxiliary  chambers  while  those  on  the  lower  right 
vent  each  end  of  the  combustion  chamber.  A  motor  case  loaded  with  propellant 
is  also  shown.  This  arrangement  provides  the  flexibility  required  to  study 


d<;R  +  oj  ) 
p  p 


2-35 


Graphite  rotor  sleeve 

(holes  in  phase  or  180°  out  of  phase) 


Figure  2-17.  Apparatus  Layout 
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Figure  2— 1 8 .  Dual  Rotating  Valve  Components 
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both  velocity  coupled  and  pressure  coupled  configurations  simply  by  changing 
the  phase  between  the  two  rows  of  holes  in  the  graphite  rotor  sleeve. 

Experimental  studies,  which  will  be  described  later  in  this  report, 
demonstrated  that  the  spacing  between  holes  and  the  relative  position  of  the 
two  rows  of  holes  in  the  rotor  are  critical.  Essentially  the  apparatus 
depends  on  "phase  locking"  the  two  valves  together  and  running  the  test  at  a 
constant  frequency.  Significant  variations  in  frequency  are  produced  if  the 
holes  in  each  row  are  not  spaced  uniformly.  The  result  is  significant  modula¬ 
tion  of  the  pressure  amplitude  and  substantial  errors  in  the  phase  angle  mea¬ 
surement.  Considerable  effort  was  required  to  define  and  correct  this  prob¬ 
lem.  This  aspect  of  the  apparatus  represented  the  major  modification  from  the 
single  rotating  valve. 

Figure  2-19  shows  a  photograph  of  the  apparatus  with  a  propellant  grain 
being  inserted  into  the  combustion  chamber.  The  steady-state  nozzle  and 
igniter  wire  inserted  in  preparation  for  a  combustion  test  are  shown  in  fig¬ 
ure  2-20.  Kistler  pressure  transducers  are  located  at  each  end  of  the  grain 
to  monitor  the  oscillating  component  of  pressure.  In  addition,  a  Taber  trans¬ 
ducer  monitors  the  average  chamber  pressure.  The  valve  at  each  end  of  the 
chamber  also  has  an  auxiliary  chamber  which  serves  as  a  phase  reference.  Each 
is  equipped  with  a  Kistler  and  Taber  transducer  to  monitor  the  oscillating 
pressure  and  the  average  chamber  pressure.  All  the  transducer  outputs  are 
recorded  with  an  FM  tape  recorder  and  played  back  through  appropriate  filters 
and  phase  meters  as  described  in  references  16  and  17.  Figure  2-21  shows  the 
instrumentation  facility  associated  with  this  apparatus. 

The  individual  slots  were  calibrated  by  using  the  clockwise-counter¬ 
clockwise  rotation  procedure  developed  for  the  single  rotating  valve.  To 
calibrate  the  slot  at  Z  =  0,  the  vent  from  the  combustion  chamber  at  Z  =  1 
was  plugged.  Pressure  coupled  calibrations  were  run  with  the  valve  rotating 
in  both  directions.  Comparisons  were  made  with  both  the  auxiliary  chambers. 
This  procedure  was  then  repeated  with  the  slot  at  Z  =  0  plugged.  Results 
from  these  two  calibrations  were  then  used  to  correct  for  phase  differences 
produced  by  tolerance  build-up. 
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Figure  2-19.  Propellant  Grain  and  Valve  Assembly 


Figure  2-20.  Assembled  Rotating  Valve  Apparatus 


2-38 


I 


I 

I 


10900-9 


Figure  2-21.  Instrumentation  Facility 
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2.3.2  Cold-Flow  Tests 

Cold-flow  tests  were  conducted  as  the  first  step  in  experimentally 
evaluating  this  apparatus.  In  these  studies,  nitrogen  was  injected  into  the 
combustion  chamber  and  the  two  auxiliary  chambers  through  individual  sonic 
chokes.  The  discharge  coefficients  of  all  three  exhaust  nozzles  were 
evaluated  by  calibration  against  a  standard  Venturi  flow  meter.  Under  these 
conditions  the  response  functions  are  zero  and  the  analysis  described  in  the 
previous  sections  can  be  used  to  predict  the  ballistics  of  all  three  chambers. 


The  first  series  of  tests  was  conducted  by  using  a  rotor  sleeve  where  the 
two  rows  of  holes  were  in  phase,  i.e..  the  pressure  coupled  configuration. 
Figure  2-22  shows  the  excellent  agreement  obtained  between  the  predicted  and 
observed  amplitudes  in  all  three  chambers.  These  tests  were  conducted  at 
frequencies  between  100  and  250  Hz  by  using  the  clockwise-counterclockwise 
method  of  reference  16. 

I 

2-39 


f 


Predicted  amplitude 


Figure  2-22.  Cold-Flow  Amplitude  Data  Pressure-Coupled  Dual 

Valve  Configuration 


129^1 


The  corresponding  phase  comparison,  shown  in  figure  2-23,  used  the  oscil¬ 
lating  pressur-  in  the  combustion  chamber  as  a  reference.  Again,  excellent 
agreement  was  obtained  between  the  predicted  and  observed  phase  differences. 
These  results  are  important  because  they  demonstrate  that  the  operation  of  the 
apparatus  is  basically  sound  when  the  two  valves  are  operated  simultaneously. 

These  initial  tests  were  conducted  by  using  only  cne  Kistler  pressure 
transducer  in  the  combustion  chamber.  The  analytical  studies  described  in  par¬ 
agraph  2.2,  however,  showed  that  pressure  measurements  at  each  end  of  the 
burner  are  required.  In  addition,  studies  at  AFRPL  showed  that  replacing 
the  graphite  slots  with  90%  Ta,  10%  W  slots  increased  their  service  life 
significantly.  Finally,  the  stainless  steel  rotor  body  under  the  graphite 
sleeve  was  replaced  with  cne  made  from  90%  Ta,  10%  W  to  increase  the  service 
life.  A  second  series  of  cold-flow  tests  was  conducted  to  insure  proper  opera¬ 
tion  with  these  modifications. 


2-UO 


Figure  2-23. 
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Cold-Flow  Phase  Data,  Pressure-Coupled  Dual  Valve  Configuration 
(Combustion  Chamber  as  Reference) 
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The  results  from  this  second  series  of  tests  are  shown  in  figures  2-2*4 
through  2-27*  They  show  excellent  agreement  between  the  predicted  and 
measured  amplitudes  and  phase  angles.  The  phase  results  show  little  differ¬ 
ence  when  the  combustion  chamber  pressure  is  measured  at  each  end  of  the 
burner.  This  substantiates  the  prediction  that  the  burner  operates  in  the 
bulk  mode  at  these  low  frequencies.  The  excellent  agreement  also  substanti¬ 
ates  the  utility  of  the  clockwise-counterclockwise  method  for  calibrating  the 
individual  slot  dimensions. 

2.3.3  Combustion  Tests 

The  initial  combustion  tests  were  conducted  to  evaluate  the  performance 
of  the  apparatus.  The  first  tests  were  conducted  with  the  standard  rotating 
valve  grains  (2.75-in.  long  by  1.0-in.  diameter  initial  port  with  a  0.25-in. 
web)  and  a  steady  state  nozzle  at  one  end.  Figure  2-28  shows  oscilloscope  pic 
tures  of  the  oscillating  pressures  at  each  end  of  the  burner  and  the  mean  cham 
ber  pressure  versus  time.  The  time  delay  between  the  rise  of  Pc  and  the 
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Blocked  Vent  at  7.  -  0  Figure  2-25.  Cold  Flow  with  Blocked  Vent  at  Z 


O  Main  Pc  at  Z  -  0  versus 
auxiliary  at  Z  =  0 


Figure  2-26.  Cold-Flow  with  Blocked  Vent  at  Z  =  1  Figure  2-27.  Cold-Flow  with  Blocked  Vent  at  z 


appearance  of  oscillations  resulted  from  delays  in  the  release  of  the  Kistler 
grounding  circuit.  Results  are  shown  for  various  filter  settings,  Q,  at  the 
fundamental  frequency  and  for  the  second  harmonic  content. 

The  most  notable  observations  are  the  large  amplitude  modulations  in  the 
oscillating  pressure.  By  imposing  the  condition  that  3e/aT  =  iXe,  the  anal¬ 
ysis  does  not  permit  evaluation  of  the  modulations.  To  investigate  their 
source,  a  series  of  tests  was  made  using  equal  steady-state  nozzles  at  each 
end  of  the  motor.  The  analysis  predicts  the  amplitude  of  the  oscillations 
should  be  extremely  snail  under  these  conditions.  The  test  results  in  figure 
2-29  show  little  change  from  the  previous  case,  which  suggests  a  serious  prob¬ 
lem  in  either  the  apparatus  or  the  analysis. 

Further  study  revealed  that  the  frequency  of  the  modulations  is  every 
twenty  cycles,  which  equals  the  number  of  holes  around  the  rotor.  Careful 
review  of  the  machining  tolerances  showed  these  rotors  had  minor  variations  in 
the  circumferential  spacing  between  holes.  New  rotors  were  machined  to  more 
exacting  tolerances .  Tests  with  these  new  rotors  showed  a  substantial  reduc¬ 
tion  in  the  magnitude  of  the  modulations,  as  shown  in  figure  2-30.  Hence,  the 
predominate  source  of  the  modulations  was  the  machining  tolerances  of  the 
graphite  rotor  sleeves. 

Additional  combustion  tests  were  then  conducted  to  obtain  quantitative 
velocity  coupled  combustion  response  function  data.  In  particular,  these  tests 
were  made  to  evaluate  the  reproducibility  of  the  data  and  the  ability  to  dis¬ 
tinguish  differences  in  propellant  formulations.  Two  nonaluminized  formula¬ 
tions  and  one  aluminized  formulation  were  selected.  The  two  nonaluminized  for¬ 
mulations,  UTP-19,933  and  19,9*12,  were  chosen  because  they  had  been  tested  in 
the  T-burner  and  rotating  valve  for  low  frequency  pressure  coupled  response. ^5 
These  data  are  shown  in  figure  2-31.  More  importantly,  the  high  frequency  com¬ 
bustion  stability  properties  had  been  studied  in  small  motor  tests.  Signifi¬ 
cant  differences  in  high  frequency  combustion  response  were  noted  between 
these  two  formulations,  which  differ  only  by  the  addition  of  0.5>  of  a  copper 
combustion  modifier.  The  burning  rates  are  nearly  identical  (±10%)  for  combus¬ 
tion  pressures  below  1,500  psi.  Since  high  frequency  combustion  stability 
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Figure  2-31.  Pressure  Response  Data  for  UTP-19,933  and  UTP-19,942 
V— 11812  13  649 

involves  gas  phase  processes,  one  could  argue  that  perhaps  these  propellants 
would  also  have  different  velocity  coupled  responses  as  well. 

Figures  2-32  through  2—34  show  the  imaginary  part  of  the  linear  velocity 
response  versus  frequency  for  UTP-19,933,  the  formulation  which  does  not  con¬ 
tain  the  additive.  Since  variations  in  propellant  burning  rate  were  noted  in 
these  tests,  figures  2-33  and  2-34  are  attempts  to  account  lor  these  varia¬ 
tions  in  the  frequency.  In  figure  2-33,  the  frequency  has  been  normalized  to 
the  burning  rate  as  suggested  by  Cohen. In  figure  2-34,  the  frequency  has 
been  normalized  to  the  square  of  the  burning  rate  as  suggested  by  thermal  wave 
combustion  models.  The  imaginary  part  of  the  response  was  selected  since  this 
part  of  the  response  determines  the  velocity  coupled  contribution  to  combus¬ 
tion  stability  in  motors. 

These  results  show  reasonable  reproducibility  from  test-to-test ,  particu¬ 
larly  when  one  considers  the  early  state  of  development  of  this  test  method 
and  the  reproducibility  difficulties  in  the  velocity  coupled  T-burner. 
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Figure  2-32.  Velocity-Coupled  Response  Function  from  Dual  Rotating 

Valve  Data  (UTP-19,933) 


Figure  2-33-  Velocity-Coupled  Response  from  Dual  Rotating  Valve  Data 
(Frequency  Corrected  by  f(r  /r),  UTP-19,933) 
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Figure  2-3^.  Velocity-Coupled  Response  Function  from  Dual  Rotating  Valve 
Data  (Frequency  Corrected  by  f(r  /r)^,  UTP-19,933) 
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Furthermore,  the  magnitude  of  the  imaginary  part  is  reasonable  when  compared 
to  previously  reported  values  obtained  in  the  T-burner.^’^  At  first  one 
might  consider  the  negative  values  to  be  suspicious  at  the  lower  frequencies, 
based  on  the  thermal  wave  analogy  to  pressure  coupled  response.  T-burner 
data  for  nonaluminized  propellants,  however,  have  also  shown  similar  behavior, 
so  these  results  are  considered  reasonable. 

Figures  2-35  through  2-37  show  the  results  for  UTP-1 9.9^2,  the  non- 
aluminized  formulation  containing  the  combustion  modifier.  Again,  burning 
rate  variations  were  noted  between  tests  and  these  effects  have  been  included 
in  the  frequency  normalization.  Again,  reasonable  reproducibility  was  found 
between  tests  and  between  different  batches  of  the  same  propellant,  particular¬ 
ly  if  one  accepts  the  frequency  normalizations  and  the  possibility  of  multiple 
peaks  in  the  response.  More  tests,  however,  particularly  at  frequencies 
between  MOO  Hz  and  600  Hz,  are  needed  to  establish  these  observations  more 
conclusively. 
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Figure  2-37.  Velocity-Coupled  Response  Function  from  Dual  Rotating  Valve 
Data  (Frequency  Corrected  by  f(r0/r)2,  UTP-19,942) 
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Comparing  the  results  for  these  two  propellants  shows  significant  differ¬ 
ences  in  the  magnitude  of  the  imaginary  part  of  the  response.  In  addition, 
there  are  important  differences  in  the  effect  of  frequency  on  the  response. 

This  suggests  the  test  method  can  distinguish  differences  in  propellant 
formulations  and  thereby  might  be  useful  for  qualitatively  ranking  propellants. 
Comparisons  between  rotating  valve  and  motor  test  results  are  needed,  however, 
to  establish  this  point  conclusively. 

Tests  were  also  conducted  with  UTP-19,360,  a  low  burning  rate  aluminized 
propellant,  since  this  class  of  formulation  is  generally  sensitive  to  erosive 
burning.  Figures  2-38  through  2-40  show  the  imaginary  part  of  the  response 
plotted  against  the  three  frequency  parameters.  These  results  are  more  diffi¬ 
cult  to  evaluate;  one  could  conclude  the  scatter  is  excessive  or  one  could  con¬ 
clude  the  response  is  characterized  by  multiple  peaks.  In  an  effort  to  clar¬ 
ify  this  point,  the  pressure  response  function  data  from  the  single  rotating 
valve  are  shown  in  figure  2-41.  If  an  analogy  between  pressure  and  velocity 
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Figure  2-39.  Velocity-Coupled  Response  Function  from  Dual  Rotating  Valve 
Data  (Frequency  Corrected  by  f(r  /r) ^UTP-19,360) 
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response  oan  be  made  from  thermal  wave  arguments,  then  one  would  conclude  the 
velocity  response  exhibits  considerable  scatter. 

Strand,  ^  however,  has  also  made  pressure  response  function  measurements 
on  this  propellant  in  the  microwave  bomb.  His  data,  which  are  shown  in  figure 
2-42,  suggest  multiple  peaks  in  the  pressure  response.  (They  also  show  a  sig¬ 
nificant  difference  with  the  rotating  valve  data.)  This  would  lead  one  to  sug¬ 
gest  the  velocity  response  should  also  show  multiple  peaks.  Thus,  contradic¬ 
tory  conclusions  can  be  drawn  by  following  this  line  of  logic  as  well.  This 
suggests  that  more  work  is  required  to  resolve  the  reproducibility  question 
for  this  aluminized  propellant. 

The  preceding  results  were  interpreted  by  using  the  linear  velocity 

response  model.  As  noted  earlier,  nonlinear  velocity  coupling  is  another  possi- 
......  .  •  •  •  -  •  *  .  .... 

ble  framework  for  interpreting  velocity  response.  This  model  would  produce  a 
high  harmonic  content  in  the  pressure  wave;  hence,  data  analyses  were  con¬ 
ducted  to  determine  the  harmonic  content  of  the  pressure  oscillations.  Ampli¬ 
tude  spectra  analysis  from  two  tests,  both  conducted  at  300  Hz  (  X  =  0 . 1 8 ) ,  are 
shown  in  figures  2-43  and  2-44.  The  data  shown  in  figure  2-43  were  obtained 
by  using  one  steady-state  nozzle,  while  the  data  in  figure  2-44  were  obtained 
by  using  identical  steady-state  nozzles  at  each  end  of  the  burner.  Additional 
frequency  spectra,  obtained  from  a  test  run  at  600  Hz,  are  shown  in  fig- 
2-45.  Several  ensembles  are  shown  simultaneously  to  demonstrate  that  the 
spectra  are  essentially  constant  during  the  entire  test. 

Several  observations  can  be  made.  First,  these  figures  show  significant 
harmonic  content.  Analysis  of  the  area  wave  form1^  shows  no  driving  of  the 
even  harmonic  by  the  valve  and  only  minimal  (-25  dB)  driving  of  the  third  har¬ 
monic.  Thus,  the  energy  at  the  even  harmonics  is  produced  either  by  the  com¬ 
bustion  process  or  by  the  fluid  mechanics  of  the  burner,  but  not  by  the 
rotating  valve  drivers. 

Calculations  were  then  made  to  determine  if  nonlinear  acoustics  or  non¬ 
linear  velocity  coupling  could  account  for  this  harmonic  content.  The  non¬ 
linear  acoustic  calculations  described  in  section  2. 2. 1.4  clearly  established 
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Figure  2-^5.  Frequency  Spectra  During  Test 


the  amplitude  of  the  second  harmonic  is  of  the  order  of  -55  dB  of  the  fundamen¬ 
tal  component.  Hence,  nonlinear  acoustics  could  not  generate  these  observed 
spectra. 

On  the  other  hand,  nonlinear  velocity  coupling  could  produce  harmonic 
amplitudes  of  -25  dB  compared  to  the  fundamental.  The  pressure  amplitudes 
generated  by  nonlinear  velocity  coupling  are  shown  in  table  2-4  for  two  values 
of  X^ .  The  values  for  Xi  =  0.1256  simulate  conditions  with  one  steady-state 
nozzle  while  the  results  for  =  0.5  simulate  identical  steady-state  nozzles 
at  each  end  of  the  burner.  These  calculations  show  the  change  in  steady-state 
nozzle  decreases  the  magnitude  of  the  second  harmonic  only  slightly.  This  is 
consistent  with  the  data  shown  as  figures  2-43  and  2-44. 

Table  2-4,  however,  does  not  explain  the  magnitude  of  the  pressure  oscil¬ 
lation  at  the  driver  frequency  shown  in  figure  2-43  and  2—^14 .  The  calcula¬ 
tions  show  the  pressure  oscillation  drops  significantly  when  X1  is  increased 
to  0.5.  Intuitively,  ere  would  expect  no  pressure  oscillations  at  the  funda¬ 
mental  under  these  conditions.  Additional  analyses  show  that  momentum  effects 
generate  the  values  shown  in  table  2-4.  The  data  shown  in  figures  2-43  and 
2-44  do  not  show  the  significant  decrease  in  amplitude  when  X-|  is  increased. 

On  the  other  hand,  the  data  shown  in  figures  2-43  and  2-45  suggest  the 
amplitude  of  the  second  harmonic  relative  to  the  amplitude  at  the  fundamental 
is  approximately  -18  to  -20  dB.  This  is  at  least  consistent  with  the  values 
shown  in  table  2-4  for  X-,  =  0.1256.  Since  table  2-4  was  calculated  by  using 
estimated  combustion  responses,  this  agreement  is  relatively  good. 


2-58 


TABLE  2-4.  PREDICTED  PRESSURE  AMPLITUDES  FOR  NONLINEAR  VELOCITY 
COUPLING  FOR  ANB  3066  WITH  LINEAR  ACOUSTICS  (Mfh  =  0.00) 

T6196 

X,  =  0.1256  X1  =  0.5 


Driver  At  Driver  At  Second  At  Driver  At  Second 


Frequency  Frequency  Harmonic  Frequency  Harmonic 


0.0598 

0.762 

X 

10~2 

0.367 

X 

10"3 

0.986 

X 

10-3 

0.432 

X 

10“ 3 

0.1197 

0.527 

X 

10"2 

0.346 

X 

1-3 

0.268 

X 

10"3 

0.353 

X 

lO-3 

0.1795 

0.452 

X 

10~2 

0.314 

X 

10-3 

0.606 

X 

10-3 

0.294 

X 

10-3 

0.2394 

0.418 

X 

10-2 

0.202 

X 

10"3 

0.783 

X 

10~3 

0.184 

X 

lO-3 

0.2992 

0.377 

X 

10~2 

0.123 

X 

10*3 

0.956 

X 

lO*3 

0.108 

X 

10" 3 

0.3590 

0.303 

X 

10"2 

0.810 

X 

10“4 

0.113 

X* 

10”-2- 

0.694 

X 

10r4 

0.4189 

0.214 

X 

10“2 

0.590 

X 

10“^ 

0.131 

X 

10-2 

0.505 

X 

10-4 
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3.0  CONCLUSIONS 


Two  analytical  models  of  the  dual  rotating  valve  were  developed  to  esti¬ 
mate  the  burner  performance  and  derive  velocity  coupled  response  functions 
from  experimental  data.  Both  models  are  based  on  one-dimensional,  linear  gas 
dynamics;  one  includes  a  linear  velocity  ooupling  model,  and  the  second  incor¬ 
porates  an  amplitude  dependent  velocity  coupling  model.  The  effects  of  pres¬ 
sure  coupling,  flow  turning,  and  particle  damping  are  also  included  in  both 
models.  Parametric  and  analytical  studies  demonstrated  that  the  equations  pre¬ 
dict  the  proper  bulk  mode  behavior  at  low  frequencies.  The  analyses  also  pre¬ 
dict  the  linear  damping  at  the  natural  acoustic  frequencies  which  agrees  with 
alternative  prediction  methods.  Finally,  excellent  agreement  was  demonstrated 
with  the  model  for  the  modulated  throat  rocket  motor. 

• 

Analytical  studies  demonstrated  the  utility  of  the  data  reduction  methods. 
For  the  linear  velocity  coupling  model,  the  data  reduction  process  was  found 
to  introduce  a  maximum  error  of  ±20%.  A  typical  error,  however,  was  found  to 
be  closer  to  ±10%.  Analyzing  the  data  with  the  nonlinear  velocity  coupling 
model  was  more  difficult.  If  the  threshold  speed  is  known  independently,  these 
methods  provide  satisfactory  results.  Separating  the  response  function  from 
the  threshold  speed,  however,  is  extremely  difficult  and  as  yet  no  satisfactory 
method  has  been  developed. 

Cold-flow  tests  confirm  the  behavior  predicted  for  the  pressure  coupled 
configuration  at  low  frequencies.  In  the  velocity  coupled  configuration,  con¬ 
trol  cf  the  phase  angle  between  the  two  valves  was  found  to  be  extremely  impor¬ 
tant.  Methods  and  specifications  were  developed  to  achieve  the  necessary  accu¬ 
racy.  Methods  for  calibrating  the  mechanical  alignment  of  the  apparatus  were 
also  demonstrated. 

Combustion  tests  were  conducted  using  both  aluminized  and  nonaluminized 
propellants.  Satisfactory  apparatus  performance  was  demonstrated  once  certain 
materials  were  improved  and  certain  design  modifications  incorporated.  Linear 
velocity  response  functions  were  derived  for  three  propellants  for  frequencies 
between  150  Hz  and  600  Hz.  Significant  differences  in  response  were  found  for 
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these  three  formulations,  thereby  demonstrating  that  the  method  can  distin¬ 
guish  propellant  formulation  changes.  Reasonable  reproducibility  was  found, 
considering  the  early  status  of  the  development  of  this  test  method.  Batch  to 
batch  reproducibility  was  also  satisfactory  at  this  early  stage  of  development 
The  frequency  of  the  oscillations  was  found  to  be  an  important  parameter.  In 
view  of  the  complexity,  nonlinear  response  functions  were  not  derived  from 
these  test  data. 

Frequency  spectra  from  the  pressure  oscillations  showed  significant  har¬ 
monic  content.  Estimates  for  various  nonlinear  contributions  showed  the 
observed  harmonic  content  was  30  dB  higher  than  the  contributions  from  non¬ 
linear  gas  dynamics.  The  predicted  nonlinear  velocity  coupling  configurations 
however,  are  similar  to  the  observed  harmonic  content.  Hence,  these  data 
strongly  suggest  the  nonlinear  velocity  coupling  model  predicts  some  of  the 
qualitative  aspects  of  the  data  which  are  observed  experimentally.  Additional 
detailed  studies  are  needed,  however,  to  provide  additional  comparisons 
before  the  validity  of  this  model  is  fully  established. 
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